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Abstract 



We consider two Large Eddy Simulation (LES) models for the approximation of large scales 
of the equations of Magnetohydrodynamics (MHD in the sequel). We study two a-models, 
which are obtained adapting to the MHD the approach by Stolz and Adams with van Cittert 
approximate deconvolution operators. First, we prove existence and uniqueness of a regular 
0^ | weak solution for a system with filtering and deconvolution in both equations. Then we study 

the behavior of solutions as the deconvolution parameter goes to infinity. The main result of 
this paper is the convergence to a solution of the filtered MHD equations. In the final section 
\ we study also the problem with filtering acting only on the velocity equation. 
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1 Introduction 



In this paper we study the equations of (double viscous) incompressible MHD 
d t u + V • (u ® u) - V • (B <g> B) + Vp - j/Au = f, 



q I d t B + V • (B <g> u) - V • (u ® B) - ^AB = 0, 

^ ' ' V-u = V-B = 0, 

u(0,x) = u (x), B(0,x)=B (x), 

where v > is the kinematic viscosity, while /i > is the magnetic diffusivity. The fields u and B 
are the velocity and the magnetic field respectively, while the scalar p is the pressure (rescaled by 
the density supposed constant here). We consider the problem in the three dimensional setting, 
and most of the technical difficulties are those known for the 3D Navier-Stokes equations (NSE). 
Examples of fluids which can be described by these equations (1.1) are for instance plasmas, liquid 
metals, and salt water or electrolytes. See Davidson [14] for an introduction to the topic. In this 
paper we aim to study the approximate deconvolution procedure (developed for turbulent flows 
by Stolz and Adams [37, 38, 1]) and especially its adaption to the MHD with the perspective of 
numerical simulations of turbulent incompressible flows, when coupled to a magnetic field. 

In the recent years, the topic of MHD attracted the interest of many researchers and, for the 
study of the question of existence, uniqueness, regularity, and estimates on the number of degrees 
of freedom, we recall the following papers [8, 10, 20, 21, 22, 23, 24, 29]. 
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Approximate Deconvolution Models (ADM) for turbulent flows without magnetic effects were 
studied in [5, 16, 27, 28]. The problem of the limiting behavior of the models when the grid mesh 
size goes to zero is already under control [16, 25, 26, 32]. On the other hand, the question of the 
limiting behavior of the solutions when the deconvolution parameter goes to infinity is a very recent 
topic, and is well-studied just for the NSE -without any coupling- in [5] (see also a short review 
in [3]). 

In the context of MHD, the topic seems not explored yet, hence we adapt here the results of [5] 
to the equations with the magnetic field and we find also some interesting unexpected variant, 
related to the applications of two different filters. Especially the equation for the magnetic field 
turns out to behave much better than that for the velocity, hence it seems not to require filtering. 

To briefly introduce the problem (the reader can find more details in the introduction of [5]), 
we recall that the main underlying idea of LES, see [4, 12, 34], is that of computing the "mean 
values" of the flow fields u = (it 1 , it 2 , it 3 ), B = (B 1 , B 2 , B 3 ), and p. In the spirit of the work started 
with Boussinesq [7] and then with Reynolds [33] , this corresponds to find a suitable computational 
decomposition 

u = u + u', B = B + B', and p = p + p', 

where the primed variables are fluctuations around the over-lined mean fields. In our context, the 
mean fields are defined by application of the inverse of a differential operator. By assuming that 
the averaging operation commutes with differential operators, one gets the filtered MHD equations 

d t u + V • (u® u) - V • (B ® B) + Vp - vAu = f , 

, d t B + V • (B®u) - V • (u® B) - uAB = 0, 

(1.2) v ' y ' _ 

Vu = V-B = 0, 

u(0, x) = uo(x), B(0, x) = B (x). 
This raises the question of the interior closure problem, that is the modeling of the tensors 

(c ® d) with eitherc, d = u or B 

in terms of the filtered variables (u, B,p). 

From this point, there are many modeling options. The basic model is the sub-grid model 
(SGM) that introduces an eddy viscosity of the form v t = Ch(x) 2 \\7u\, which may be deduced 
from Kolmogorov similarity theory (see [12]), where h{x) denotes the local size of a computational 
grid, and C is a constant to be fixed from experiments. This model, that already appears in Prandtl 
work [31] with the mixing length I instead of h(x), was firstly used by Smagorinsky for numerical 
simulations [36] . This is a very good model, but introduces numerical instabilities in high-gradient 
regions, depending on the numerical scheme and potential CFL constraints. 

Among all procedures to stabilize the SGM, the most popular was suggested by Bardina et 
al. [2], which reveals being a little bit too diffusive and underestimates some of the resolved scales, 
that are called "Sub Filter Scales" (SFS) (see for instance [11, 19]). Then the model needs to 
be "deconvolved" to reconstruct accurately the SFS. Hence, many options occur here, too. In 
the present paper we study the Approximate Deconvolution Model (ADM), introduced by Adams 
and Stolz [37, 1], who have successfully transferred image modeling procedures [6] to turbulence 
modeling. 

From a simplified and naive mathematical viewpoint, this model, which uses similarity proper- 
ties of turbulence, is defined by approximating the filtered bi-linear terms as follows: 

(c® d) - (D N (c) ® D N (d) . 

Here the filtering operators Gi are defined thanks to the Helmholtz filter (cf. (2.1)-(2.2) below) by 
Gi(u) = u, G2(B) = B, where Gi := (I — afA)^ 1 , i — 1,2. Observe that we can then have two 
different filters corresponding to the equation for the velocity and for that of the magnetic field. 
There are two interesting values for the couple of parameters (a>i, 02) G K + x R + : 

1. a\ = ot2 > 0. In this case the approximate equations conserve Alfven waves, see [22]; 

2. «i > 0, a2 = 0, which means no filtering in the equation for B. 
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The deconvolution operators Djv 4 are defined through the van Cittert algorithm (2.10) and the 
initial value problem that we consider in the space periodic setting is: 

d t vr + V • d (D Nl (w) ® D Nl (w)) - V • d (£>at 2 (b) ® D N2 (b)) + Vg - i/Aw = GJ, 
d t b + V- G 2 (Djv, (w) ® D W2 (b)) - V • G 2 (Dat 2 (b) ® (w)) - ^Ab = 0, 
(1.3) Vw = Vb = 0, 

w(0, x) - G lUo (x), b(0, x) = G 2 B (x), 
a\ > 0, a 2 > 0. 

As usual, we observe that the equations (1.3) are not the equations (1.2) satisfied by (u, B), but 
we are aimed at considering (1.3) as an approximation of (1.2), hence w ~ G±u and b ~ G 2 B. 
This is mathematically sound since, at least formally, 

D Ni -> Ai := I - a 2 A in the limit Ni — > +oo, 

hence, again formally, (1.3) will become the filtered MHD equations (1.2). The existence and 
uniqueness issues have been also treated (even if without looking for estimates independent of Ni) 
in [22, 21] (for arbitrary deconvolution orders). What seems more challenging is to understand 
whether this convergence property is true or not, namely to show that as the approximation 
parameters Ni grow, then (as recently proved for the Navier-Stokes equations in [5]) 

w — > Giu, b — > G 2 B, and q — > G\q. 

We prove that the model (1.3) converges, in some sense, to the averaged MHD equations (1.2), 
when the typical scales of filtration ai remain fixed. Before analyzing such a convergence, we need 
to prove more precise existence results. To this end we follow the same approach from [5], which 
revisits the approach in [16] for the Navier-Stokes equations. To be more precise, the main result 
deals with a± > and a 2 > 0. We first prove (cf. Thm. 3.1) existence and uniqueness of solutions 
(wjv,bs,gAr) of (1.3), with N = (Ni,N 2 ), such that 

WJV , h N eL 2 ([0,T]; J ff 2 (T3) 3 )ni oo ([0,T];if 1 (T 3 ) 3 ), 
q N G L 2 ([0, T]; W^ 2 (T 3 )) R L 5/3 ([0, T]; W 2 ^ 3 (T 3 )), 

and our main result is the following one. 

Theorem 1.1. Let a\ > and a 2 > 0; then, from the sequence {(wjv, bjv, ^a/OIa^n 2 ; one can 
extract a (diagonal) sub-sequence (still denoted {(wjv, bjv, <7at)}at(en 2 J 

wat^w j weakly in L 2 ([0,T}; H 2 (T 3 f) 

h N ^b [ weakly* in L°°([0, T]; H\T 3 ) 3 ), 



strongly in L p ([0, T]; i? 1 (T 3 ) 3 ), VI <p < +oo, 



W^r — > W 

b a^ ^ b 

q N -> q weakly m L 2 ([0, T}; W 1,2 (T 3 )) n L 5 / 3 ([0, T]; VK 2 ' 5 / 3 (T 3 )), 



such that the system 

<9 t w + V • Gi(AiW (g> Ajw) - V • Gi(A 2 b ® A 2 h) + Vq - vAw = Gif, 

Vw = Vb = 0, 

(1.4) 

d t h + V • G 2 (A 2 b <g> A x w) - V • G 2 (Aiw ® A 2 b) - ^Ab = 0, 
w(0, x) - G lUo (x), b(0, x) = G 2 B (x), 

holds in the distributional sense. Moreover, the following energy inequality holds: 
(1.5) i^(||A 1 w|| 2 + ||A 2 b|| 2 )+ i ,||V^ 1 w|| 2 + M ||VA 2 b|| 2 < (f,A x w). 
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As a consequence of Theorem 1.1, we deduce that the field (u, B,p) = (Aiw, A%b, A%q) is a 
dissipative (of Leray-Hopf 's type) solution to the MHD Equations (1.1). 

Remark 1.1. Following the work carried out in [17] about ADM without coupling, we conjecture 
that the error modeling in the case of the approximate deconvolution MHD is of order _/V~ 1//2 . 

Remark 1.2. The question of the boundary conditions is the great challenge in LES modeling, 
see [35] for a general review. This is why, either from physical or mathematical viewpoint, the- 
oretical issues are raised in the case of periodic boundary conditions, although the reality of such 
boundary conditions may be controversial. In a paragraph in [38], the authors outline a possible 
numerical discrete algorithm for a deconvolution procedure, by the finite difference method. This 
is proposed in the case of homogeneous boundary conditions, but there is no mathematical analysis 
about this method, which still remains an open problem. 



Plan of the paper. In Sec. 2 we introduce the notation and the filtering operations. In Sec. 3-4 
we consider the model with the double filtering with non-vanishing parameters OLi and then we 
study the limiting behavior as Ni — > +oo. Finally, in Section 5, we treat the same problems in 
the case a.\ > and cxi = 0. Since most of the calculations are in the same spirit of those in [5], 
instead of proofs at full length we just point out the changes needed to adapt the proof valid for 
the NSE to the MHD equations. 



2 Notation and Filter/Deconvolution operators 

This section is devoted to the definition of the functional setting which we will use, and to the 
definition of the filter through the Hclmholtz equation, with the related deconvolution operator. 
All the results are well-known and we refer to [5, 27, 28] for further details. We will use the 
customary Lebesgue LP and Sobolev W k,p and W 11,2 = H s spaces, in the periodic setting. Hence, 
we use Fourier series on the 3D torus T3. Let be given L G := {x G R : x > 0}, and define 
SI :=]0, L[ 3 C R 3 . We denote by (ei, e2, 63) the orthonormal basis of R 3 , and by x := (x\, X2,x 3 ) G 
R 3 the standard point in R 3 . We put % := 2ttS 3 / 'L and T 3 is the torus defined by T 3 := (R 3 /T 3 ). 
We use || • || to denote the L 2 (Ta)-norm and associated operator norms. We always impose the zero 
mean condition on the fields that we consider and we define, for a general exponent s > 0, 

H s := \ w : T 3 -4 R 3 , w G H S (T 3 ) 3 , V • w = 0, / wdx = 



For w G H s , we can expand the fields as w(x) = XkeT* w ke +lk x , where k G T3 is the wave- 
number, and the Fourier coefficients are Wk := pl^- f Ta w(x)e _ik ' x (ix. The magnitude of k is defined 
by k := |k| = {|fci| 2 + |fc 2 | 2 + |fc 3 | 2 }'- We define the H s norms by ||w|| 2 := Eker,* |k| 2s |w k | 2 , 
where of course ||w||q := ||w|| 2 . The inner products associated to these norms are (w, v)h, :== 
EkeT 3 * |k| 2s ^k • v k , where v k denotes the complex conjugate of Vk- To have real valued vector 

fields, we impose w_k = Wk for any k G T 3 * and for any field denoted by w. It can be shown (see 
e.g. [15]) that when s is an integer, ||w|| 2 :— ||V s w|| 2 and also, for general s G M, (H s )' = H_ s . 

We now recall the main properties of the Helmholtz filter. In the sequel, a > denotes a given 
fixed number and for w G H s the field w is the solution of the Stokes-like problem: 



(2.1) 



—a 2 Aw + w + V7r = w in T3, 
V • w = in T3 , 



7r dx = 0. 

T 3 



For w G H s this problem has a unique solution (w, it) G H s+ 2 x H S+1 (T 3 ), whose velocity is 
denoted also by w = G(w). Observe that, with a common abuse of notation, for a scalar function 
X we still denote (this is a standard notation) by x the solution of the pure Hclmholtz problem 

(2.2) Ax:= ~a 2 Ax + x = X in T 3 . 
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In particular, in the LES model (1.3) and in the filtered equations (1.2)— (1.4), the symbol " " 
denotes the pure Helmholtz filter, applied component-wise to the various tensor fields. 
We recall now a definition that we will use several times in the sequel. 

Definition 2.1. Let K be an operator acting on H s . Assume that e~ jkx are eigen-vectors of K 
with corresponding eigenvalues K^. Then we shall say that is the symbol of K . 

The deconvolution operator Dm is constructed thanks to the Van-Cittert algorithm by Dn '■= 
y^_ (I — G) n . Starting from this formula, we can express the deconvolution operator in terms of 
Fourier series Djy(w) = X^ker* ^Jv(k)wke +4kx , where 

<") 6 »« = E(Tfiw) =(1+0,HV " [Trim ■ 

The basic properties satisfied by Dn that we will need are summarized in the following lemma. 

Lemma 2.1. For each N £ N the operator Dn ■ Hs — > H s is self-adjoint, it commutes with 
differentiation, and the following properties hold true: 



(2.4) 1 < D N (k) < N + 1 Vke7T ? ; 

(2.5) D N (k) « (N + l) ^ff 2 for large |k|; 

(2.6) lim Djv(k) = N + 1 for fixed a > 0; 

|k|-s-+oo 

(2.7) D N (k) < (l + a 2 |k| 2 ) Vke7i, a > 0; 

(2.8) the map w i— > Djv(w) is an isomorphism s.t. ||-Djv||h s = 0(./V + 1) Vs > 0; 

(2.9) lim D w (w) = iw m H s Vs 6 R and w £ H s+2 . 

TV— > +00 



All these claims follow from direct inspection of the formula (2.3) and, in the sequel, we will 
also use the natural notations Gi :— A^ 1 := (I — afA)^ 1 and 

Ni 

(2.10) D Ni :=Y^(l-Gi) n , i = l,2. 

n—0 

3 Existence results 

In order to be self-contained, we start by considering the initial value problem for the model (1.3). 
In this section, N\, N2 £ N are fixed as well as ct\ > 0, 0/2 > 0, and we assume that the data are 
such that 

(3.1) u ,B eH , f eL 2 ([0,T] xT 3 ), 

which naturally yields G1U0, G2B0 G H2, G\i € L 2 ([0, T]; H 2 ). We start by defining the notion of 
what we call a "regular weak" solution to this system. 

Definition 3.1 ("Regular weak" solution). We say that the triple (w, b, q) is a "regular weak" 
solution to system (1.3) if and only if the three following items are satisfied: 
1) Regularity: 



(3.2) w, b e L 2 ([0, T]; H 2 ) D C([0, T]; Hi), 

(3.3) d t w,d t beL 2 ([0,T];B. o ) 

(3.4) qeL 2 ([0,T];H x (T 3 )), 

2) Initial data: 

(3.5) lim||w(i,-)-GiUo|| Hl = 0, lim ||b(i, ■) - G 2 B || Hl = 0, 
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3) Weak Formulation: For all v,h G L 2 ([0, T]; iJ 1 ^) 3 ) 



(3.6) 



[ d t w-v- [ I G x (D Nl {w)®D Nl {vf)):Vv 
Jt 3 Jo Jt 3 

+ f f G 1 (D N2 (b)®D N2 {b)):Vv+ [ [ Vq 
Jo Jt 3 Jo Jt 3 

+ v [ f Vw:Vv= / / Gi(f)-v, 
Jo Jt 3 Jo Jt 3 



(3.7) 



/ / d t b h- f ( G 2 (%(b)®%(w)):Vh 
Jo Jt 3 Jo Ji 3 

+ [ [ G 2 (D m (w) ® D N2 (b)) : Vh + /x / f Vb : Vh = 0. 

Jo Jt 3 Jo Jt 3 



Observe that, for simplicity, we suppressed all c?x and dt from the space-time integrals. With 
the same observations as in [5], one can easily check that all integrals involving -DjviW and Djy 2 b 
in (3.6)-(3.7) are finite under the regularity in (3.2)-(3.3). We now prove the following theorem, 
which is an adaption of the existence theorem in [5] and at the same time a slightly more precise 
form of the various existence theorems available in literature for doubly viscous MHD systems. 

Theorem 3.1. Assume that (3.1) holds, < a; e R and ]Vj € N, i = 1,2, are given and fixed. 
Then, problem (1.3) has a unique regular weak solution. 

In the proof we use the usual Galerkin method (see for instance the basics for incompressible 
fluids in [30]) with divergence- free finite dimensional approximate velocities and magnetic fields. 
We also point out that Theorem 3.1 greatly improves the corresponding existence result in [22] and 
it is not a simple restatement of those results. Some of the main original contributions are here 
the estimates, uniform in N, that will allow later on to pass to the limit when Ni —5- +oo. 

Proof of Theorem 3.1. Let be given m G N* and define V m to be the following space of real valued 
trigonometric polynomial vector fields 

V m :={w£Hi: f w(x) e~ ik - x = 0, Vk with |k| > to}. 

Jt 3 

In order to use classical tools for systems of ordinary differential equations, we approximate the 
external force f with f\/ m by means of Friederichs mollifiers. Thanks to the Cauchy-Lipschitz 
Theorem, we can prove existence of T m > and of unique C 1 solutions w m (i,x) and b m (i,x) 
(belonging to V m for all t G [0, T m [) to 



(3.8) 



(3.9) 



/ <9 t w m • v - / Gi(%(w m ) ® D Nl (w m )) : Vv 
Jt 3 Jt 3 

+ I Gi(%(b m ) ® D N2 (b m )) : Vv 
Jt 3 

+ is Vw m : Vv = / Gi(fi/ m ) • v, 

Jt 3 Jt 3 

/ d t b m h G 2 {D N2 (b m ) <g> D Nl (w m )) : Vh 
Jt 3 Jt 3 

+ G 2 (D Nl (-w m ) ® D N2 (b m )) : Vh + /i / Vb m : Vh = , 

Jt 3 Jt 3 



for all v, h G L 2 ([0, T]:V m ). 

Remark 3.1. Instead of (w m ,b m ), a more precise and appropriate notation for the solution of 
the Galerkin system would be {v/ m ,N 1 ,N2,a 1 ,a 2 i^ ) m,N 1 ,N2,a 1 ,a 2 )- We are asking for a simplification, 
since in this section Ni and on are fixed and the only relevant parameter is the Galerkin one m G N*. 



G 



The natural and correct test functions to get a priori estimates are AiDjVi(w m ) for the first 
equation and A 2 DN 2 (h m ) for the second one. Arguing as in [5], it is easily checked that both are 
in V m . Since Ai,A 2 are self-adjoint and commute with differential operators, it holds: 



Moreover, 



Gi(%(w m )®%(w m )) : V(A x D Nl (vi m ))<ht = 0, 
/ G 2 {D N2 (b m ) <g> D N2 (b m )) : V(A 2 D N2 (b m )) dx = 0. 

/ G\ (Dn 2 (b m ) ® Dn 2 (b m )) : V(AiD Nl (w m )) dx 
Jt 3 

- / G 2 (-Div 2 (b m )(g)£> iVl (w ro )) : V(A 2 £>iv 2 (b ro ))dx 
+ / G 2 (D Nl (w m ) ® £>jv,(b m )) : V(A 2 Djv 2 (b m )) dx 

= - / (£>jv 2 (b ro ) • V)Djv 2 (b TO ) • %(w m ) dx 

JT 3 

+ / {D Nl (w m ) • V)-Dat 2 (b m ) • D N2 (b m ) dx 

Jt 3 

- / (D N2 (b m ) ■ W)D Nl (w m ) • Djv 2 (b m ) dx = . 
Jr 3 

Summing up the equations satisfied by w m and b m , using standard integration by parts and 
Poincare's inequality combined with Young's inequality, we obtain 



AlDJr(w m )(t, 



(3.10) 



\AiD% 2 (b m )(t, 
ft 



\VAlDj fi (w m )\\ 2 dT+ / \\VAlDl 2 {h m )\\ 2 dT 
o Jo 

<C(||u ||, ||Bo||, i/- 1 ||f|| La([0 ,ri.H_ 1 )). 



which shows that the natural quantities under control are yl[D^(w m ) and A|Z)^ 9 (b m ). 

Since we need to prove many a priori estimates, for the reader's convenience we organize the 
results in tables as (3.11). In the first column we have labeled the estimates, while the second 
column specifies the variable under concern. The third one explains the bound in terms of function 
spaces: The symbol of a space means that the considered sequence is bounded in such a space. 
Finally, the fourth column states the order in terms of a, m and N for each bound. 



(3.11) 



Label 


Variable 


bound 


order 


a) 


A\D\ (w m ), AlD^{h m ) 


L°°([0,T] 


H )ni 2 ([0,T] 


Hi) 


0(1) 


b) 


< 2 (w m ), D^(b m ) 


L°°([0,T] 


H o )nL 2 ([0,T] 


Hi) 


0(1) 


c) 


< 2 (w m ), <(b m ) 


L°°([0,T] 


H!)nL 2 ([o,r] 


H 2 ) 


0(0 


d) 


w m, b m 


L°°([0,T] 


H )ni 2 ([0,T] 


Hi) 


0(1) 


e) 


w m, b m 


L°°([0,T] 


H x )nI?([0,T] 


H 2 ) 


0(a- L ) 


/) 


DjVi(w m ), Dn-2 (b m ) 


L°°([0,T] 


H )nL 2 ([0,T] 


Hi) 


0(1) 


9) 


%(w m ), D N2 (h m ) 


L°°([0,T] 


H!)nL 2 ([o,r] 


H 2 ) 


0(^E) 


h) 


9 t w m , d t h m 


£ 2 ([0,T];H ) 


0(0. 



In the previous table, a = o?i for w m , a — a 2 for b m , while in /i) we can take a := min{ai,Q! 2 } 
for both w m and b m . 

Proof of (3.11-a) - This estimate follows directly from (3.10). Notice also that since the 
operator A i 2 D^. has for symbol p^. k < 1, then \\AID^ Gifi/ m || < C||f|| and also 



|A?£>| r .P m G i a|| = 



A?DlGiB\\ < \\A?DIM\ < H, 
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which will be used with a = Uq, Bo. 

Proof of (3.11-b)-(3.11-c) — Let v <E H 2 . Then, with obvious notations, one has 

Il4v|| 2 = ]T(l + a 2 |k| 2 )|v k | 2 = ||v|| 2 + a 2 ||Vv|| 2 . 
ker 3 * 

11 11 
It suffices to apply this identity to v = D^(w m ), £»^(b m ) and to v = a j CJ i (w m ), c^D^ (b m ) 

(i = 1, 2, 3) in (3.10) to get the claimed result. 

Proof of (3.11-d)-(3.11-e)-(3.11-f) —These are direct consequence of (3.11-a)-(3.11-b)-(3.11-c) 
combined with (2.4). 

Proof of (3.11-g) — This follows directly from (3.11-e), together with (2.4). 

Remark 3.2. One crucial point is that (3.11-g) is valid for each N — (iVj, N2), but the bound 
may grow with Ni. 

Proof of (3.11-h) -Let us take 9 t w m ,<9 f b m £ V m as test vector fields in (3.8). We get 

||d t w m || 2 + / Gi(V • [D Nl (w m ) (8 D Nl (yr m )]) ■ d t w m 
Jt 3 

- I Gi(V- [D Na (b m )®D Na (b m )]) •5cw m + ~||Vw m || 2 = / G 1 i 1/m -d t w m , 
Jt 3 1 °X Jj 3 

||5 t b m || 2 + f G 2 (V-[D N2 (b m )^D Nl (w m )])-d t b m 

- / G 2 (V-[D Nl (w m )®D N2 (b m )})-d t b m + ^\\Vb m \\ 2 = 0. 

To estimate the time derivative, we need bounds on the bi-linear terms 

Ajv, m := Gi V • (D Nl (w m ) ® Djvj (w m )) , 
Bat,™ := Gi V • (D N2 (b m ) ® (b m )) , 
Cat,™ := G 2 V • (fl^fw™) ® £>;v 2 (b m )). 

Even if we have two additional terms, this can be easily done as in [5] by observing that, by 
interpolation inequalities, both DN 1 (w m ) and D^ 2 (b m ) belong to L 4 ([0, T]; L 3 (T3) 3 ). Therefore, 
by observing that the operator (V-) o Gi has symbol corresponding to the inverse of one space 
derivative, it easily follows that Ajv im ,Bjv,m, Cat,™ £ £ 2 ([0,T] x T 3 ) 3 . Moreover, the bound is of 
order 0(a~ 1 ) as well. 

From the bounds proved in (3.11) and classical Aubin-Lions compactness tools, we can extract 
sub-sequences {w m , b m } mS N converging to w, b £ L°°([0, T]; Hi) n L 2 ([0, T]; H 2 ) and such that 

weakly in L 2 ([0, T}; H 2 ), 

strongly in L"([0,T];Hi), VpG [l,oo[, 

weakly in L 2 ([0, Tj; H ). 

This already implies that (w, b) satisfies (3.2)-(3.3). From (3.13) and the continuity of Dn ( in 
H s , we get strong convergence of DN 1 (w m ), D^ 2 (b m ) in L 4 ([0,T] x T 3 ), hence the convergence of 



(3.12) 
(3.13) 
(3.14) 



w m w 
b m -> b 

w m ->• w 
b m ->• b 
<9 t w m ->• <9 t w 
<9 t b m ->• <9 t b 
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the corresponding bi-linear products in L 2 ([0, T] x T3). This proves that for all v, h £ L 2 ([0,T];Hi) 



(3.15) 



(3.16) 



Jt 



[ 9 t w v- / / <rVi(/) A| (wj ./) Al (wii : Vv 
11 Jt 3 

r 

+ 



I f G!(%(b)®%(b)):Vv 

JO JT 3 

+ u [ / Vw:Vv= f [ Gi(f)-v, 
Jo Jt 3 Jo Jt 3 



f d t b h f f G 2 (D N2 (b)®D Nl {w)):Vh 

:i JT 3 JO JT 3 



+ f f G 2 (D Nl (w) ®D N2 {b)) : Vh + /i f f Vb : Vh 

Jo Jt 3 Jo Jt 3 



= 0. 



To introduce the pressure, observe that taking the divergence of the equation for w, we get 



(3.17) 



Aq = V -GJ+V - An, 



for An ■— — Gi [V- (£>AT 1 (w)(g)I?jv 1 (w)) — V- (Z?jv 2 (b)®-DAf 2 (' :) ))] ■ A fairly standard application of 
De Rham's Theorem shows existence of q, and the regularity of An yields q 6 L 2 ([0, T]; H 1 ^^)). 

The meaning in which the initial data are taken is completely standard and we end the proof 
by showing uniqueness: Let (wi, bi) and (W2, be two solutions corresponding to the same data 
(u ,B ,f) and let us define, as usual, W :— wi — w 2 and B := bi — b 2 . By standard calculations 
(mimicking those employed in [5]), we get 

i|(ii44 i (w)ii 2 + ii44 2 (^)iiVHiv44 1 (w)ii 2 +Miiv44 2 (^)ii 2 

= / {D N2 {B)-W)D N2 {b x )-D Nl {W)- [ (£ Wl (W)-V)Av 1 (w 1 )-.D J v 1 (W) 
jt 3 Jt 3 

+ [ (Dn 2 (B) • V)-Djv x (wi) • Dn 2 (B) — f (D Nl (W) •VjDjv.fbi) • D Na (B), 
JT 3 Jt 3 

< 2\\D N2 (B)\\ L 4D Nl W\\ L 4^D N2 (b 1 )\\ L , + HAvxWII^IlVAviCwi)!!^, 

+ \\D N2 (B)\\ 2 L *\\VD Nl (-w 1 )\\ L 2 

<2||I>* a (B)||V*||iJ Wl ( W )||V4||v^^^ 

+ ||^ 1 (W)|| 1 / 2 ||V^ 1 (W)|| 3 / 2 ||V^ 1 (w 1 )|| + ||Z» J v a (B)|| 1 /2|]vX» J v 3 (iB)|| 3 / 2 ||V J D J v 1 (w 1 )|| . 

By using H-DjvJ = (iVj + 1), the bound of Wi,bi in L°°([0, T]; Hi), and Young's inequality, we 
obtain 



ld_ 

2dt 



(\\AID-^(W)\\ 2 + \\AlDl 2 {B)\\ 2 ) + ^\\VA!D Ni (W)\\ 2 + %\\VA2D* a (B)f 
< C(N X + l) 4 (sup HVwill 4 ) [-1 \\AlD Ni (W)f + 4l44 2 (£)f 

+c(n 2 + i) 4 (sup nvbiii 4 ) —L—fli^D^o*)!! 2 + ||4dJ 2 (s)|| 2 ]. 



In particular, we get 

' ' ^\AfD Ni (W)f + \\a\d1 2 {BW) < M[\\AlD Ni (W)\\ 2 + ||4^ a (B)|| a ] 



2dt 



where 



M := C(max{-, -}f \{N\ + l) 4 f sup ||V Wl || 4 ) + (JV a + l) 4 f sup ||Vb x 



1 1 

V \l~ ' L \ (> ' \ t >0 

Since the initial values W(0) = B(Q) are vanishing, we deduce from Gronwall's Lemma that 
4-Djvj (W) = 4- D jv 2 (#) = and wc conclude that W = B = 



□ 
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Remark 3.3. The same calculations show also that the following energy equality is satisfied 

i|(ii44 i (w)ii 2 + ii44jb)f)+Hiv44 i (w)ii 2 + A1 ||v44 2 (b)ii 2 
= (4^l 1 (Gif),4^l 1 (w)). 

As we shall see in the sequel, it seems that it is not possible to pass to the limit N — > +00 directly 
in this "energy equality" and some work to obtain an "energy inequality" is needed. 

4 Passing to the limit when N — > oo 

The aim of this section is the proof of the main result of the paper. For a given N g N, we denote by 
(wat, bjv, qiy) the unique "regular weak" solution to Problem 1.3, where N = min{iVi, N 2 } — > +00. 
For the sake of completeness and to avoid possible confusion between the Galerkin index m and 
the deconvolution index N, we write again the system: 

d t w N + V ■ d{D Nl (wjv) ® D Nl (w N )) - V ■ d (D N2 (b N ) ® D N2 (b N )) 

+ X7q N - vAw N = GJ in [0,T]xT 3 , 

d t h N + V ■ G 2 (D N2 (b N ) <E) D Nl (w N )) - V ■ G 2 (D Nl (w N ) <E> D N2 (b N )) 

- fiAb N = in [0,T] x T 3 , 

V-wat = V-bAr = in [0,T]xT 3 , 
(w JV ,b iV )(0,x) = (GiUo ) G 2 Bo)(x) in T 3 . 

More precisely, for all fixed scales a%, a 2 > 0, we set 

wtv = lim w mi jvj ,Ar 2 , ai ,a 2 
m— f +00 

and similarly for b^r. 

Proof of Thm. 1.1. We look for additional estimates, uniform in N, to get compactness properties 
about the sequences {Da^ (w^v), Djv 2 (b.iv)}./veN an d {w^v, bjvjwGN- We then prove strong enough 
convergence results in order to pass to the limit in the equation (4.1), especially in the nonlinear 
terms. With the same notation of the previous section, we quote in the following table the estimates 
that we will use for passing to the limit. The Table (4.2) is organized as (3.11) and a — min{ai, a 2 }- 



Label 


Variable 


bound 


order 


a 


wat, bjv 


L oo ([0,T];H )ni 2 ([0,T];H 1 ) 


0(1) 


b 


wat, b N 


L oo ([0,T];H 1 )nL ;a ([0,T];H 2 ) 


0(a- L ) 


c 


D Nl (-w N ), D N2 (b N ) 


L oo ([0,T];H )n J L 2 ([0,T];H 1 ) 


0(1) 


d 


9 t WAr, d t b N 


L 2 ([0,T] x T 3 ) 3 


0(a- L ) 


e 


qN 


L 2 ([0,T];i/ 1 (T 3 ))nL 5 / a ([0,T];l^ 2 ' 5 / 3 (T3)) 


0{a- r ) 


f 


d t D Nl (w N ), d t D N2 (b N ) 


i 4/3 ([0, T]; H_i) 


0(1) 



Estimates (4.2-a), (4.2-b), (4.2-c), and (4.2-d) have already been obtained in the previous section. 
Therefore, we just have to check (4.2-e) and (4.2-f). 

Proof of (4.2-e) — To obtain further regularity properties of the pressure we use again (3.17). 
We already know from the estimates proved in the previous section that An & L 2 ([0,T] x T3) 3 . 
Moreover, classical interpolation inequalities combined with (4.2-c) yield Dm 1 (wjv), Djy 2 (b^) € 
£10/3([ Q)T ] x T3 )_ Therefore, A N € L 5 / 3 ([0, T]; W /1,5 / 3 (T 3 )). Consequently, we obtain the claimed 
bound on q^- 
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Proof of (4.2-f) — Let be given v, h £ £ 4 ([0, T]; Hi). We use Dn 1 (v), D^ 2 (h) as test functions. 
By using that 3fW, d t b £ L 2 ([0,T] x T3) 3 , Dpj i commute with differential operators, Gi and Dj^ i 
are self-adjoint, and classical integrations by parts, we get 

(dtwjv, £>jvi(v)) = (d t D Nl (-w N ),v) 

= v(Aw N , D Nl (v)) + (D Nl (wjv) ® D Nl (wjy), G x D Nl (Vv)) 
- (D N2 (b N ) ® £> Wa (bjv),GiDjVx(Vv)) + (D^CdfJ.v), 



(c^, £>w a (h)) = (atl>jv 2 (bjv),h) 

= MAbjv, £>at 2 (h)) + (L>jv 2 (bjy) ® (wjy), G 2 £>jv 2 (Vh)) 
- (D Nl (w N ) ® Djv 2 (bjv), G 2 £»jv 2 (Vh)) . 

We first observe that 

KAwjv.^Cv))! = |(V£) J v(w J v),Vv)| < a x (t)||vl|i, 
\(Ab N ,D N2 (h))\ = \(VD N2 (b N ),Vh)\ < Ci(t)||h||i, 

and that the L 2 ([0, T]; i/ 1 (T 3 ) 3 ) bound for D Nl (w N ), D N2 (b N ) implies that d(t) £ L 2 ([0,T]), 
uniformly with respect to N £ N. Therefore, when we combine the latter estimates with the 
properties of we get, uniformly in N, 

\(d t D Nl (w N ),v)\ + \(d t D N2 (b N ),y)\ 

< (vCi(t) + C a (*))||v||i + (/iCx(i) + C a (t))||h||i + ||f(i, -)|| G L 4/3 (0,T). 

From the estimates (4.2) and classical rules of functional analysis, we can infer that there exist 

w,bGL oo ([0,r];H 1 )nL 2 ([0,T];H 2 ), 

Zl ,z 2 eL oo ([0,T];H )n J L 2 ([0,T];H 1 ), 

q £ L 2 ([0,T}; H\T 3 )) D L 5 / 3 {[0,T};W 2 ' 5/3 (T 3 )) 

such that, up to sub-sequences, 



(4.3) 



wjy — > w 
b N -> h 



9 4 wat — > <9 t w 
d t h N ->■ d t b 

D Nl (w N ) -> zi 

£>AT 2 (bAr) -> Z 2 

d t D Nl (w N ) —■ d t zx 
d t D N2 (b N ) -» d t z 2 



qN -> q 

We notice that 

D Nl (-w N ) <g> D Nl (vr N ) 
(4.4) £»jv 2 (b w ) ® L»Ar 2 (b w ) - 

DnA^n) ® £»7v 2 (bAr) 



weakly in L 2 ([0, T]; H 2 ), 
weakly* in L°°([0, T]; H x ), 
strongly in £P([0,T];Hi) 

weakly in L 2 ([0, T] x T 3 ), 

weakly in L 2 ([0, T]; Hi), 
weakly* in L°°([0, T]; H ), 
strongly in L p ([0,T] x T 3 ) 3 



V> G [l,oo[, 



v P e[i,io/3[, 



weakly in L 4/3 ([0, T]; H_i), 

weakly in L 2 ( [0, T] ; iJ 1 (T 3 ) ) n i 5/3 ( [0, T] ; W 2 > 5 / 3 (T 3 ) ) . 

-^zi®zi strongly in L p ([0,T] x T 3 ) 9 V>G[l,5/3[ 
->• z 2 ® z 2 strongly in L*([0, T] x T 3 ) 9 Vp G [1, 5/3[, 
^zi®z 2 strongly in L p ([0,T] x T 3 ) 9 V>e[l,5/3[, 



while all other terms in the equation pass easily to the limit as well. By using the same identification 
of the limit used in [5], we can easily check that zi = A x w and z 2 = A 2 b, ending the proof. □ 
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By using well established results on semicontinuity and adapting calculations well-known for 
the NSE, we can prove that the solution (w, b) satisfies an "energy inequality." 

Proposition 4.1. Let be given u ,B G Ho, f G L 2 ([0, T]; H ), and let {(wjy, bjy, <Z;v)}iV6N be 
a (possibly relabelled) sequence of regular weak solutions converging to a weak solution (w, b, q) 
of the filtered MHD equations. Then (w, b) satisfies the energy inequality (1.5) in the sense of 
distributions (see also [13, 18, 39]). This implies that (w,b) is the average of a weak (in the 
sense of Leray-Hopf) or dissipative solution (u, B) of the MHD equation (1.1). In fact, the energy 
inequality can also be read as 

i|(||u|| 2 + ||B|| 2 ) + ^||Vu|| 2 + M ||VB|| 2 < (f, u). 

Proof. The proof is a straightforward adaption of the one in [5] . We start from the energy equality 
for the approximate model as in Remark 3.3 and we observe that the same arguments as before 
show also that 



Z?] v / i 2 (w Ar )^< 2 (w)~ 



weakly in L 2 ([0,T];Hi 



Next, due to the assumptions on f, we have A 1 1 ^ 2 D^ 2 f — > f strongly in L 2 ([0,T];Ho) and, since 
for all N e N we have wjv(O) = Gi.u(0) G H 2 and bAr(O) = G 2 b(0) G H 2 , we get 

l(\\Al /2 D^(w N )(0)f + ||4 /2 ^ 2 2 (b A ,)(0)|| 2 ) + jf* (A- 1/2 < 2 (f),4 /2 < 2 (w w )) ds 

[ (f,A-Lw)ds. 
Jo 



^»l(||A lW (0)|| 2 + [[A 2 b(Q)|[ a ) 



Next, we use the elementary inequalities for liminf and limsup to infer that 
1 



Hmsup-dlA^^^w^)^)!! 2 + \\Al /2 D^(b N )(t)\\ 2 ) 
+ hminf (uj* \\WAl /2 D]/ 2 (w N )(s)\\ 2 ds + n^ \\VAl /2 D l £(b N )(s)\\ 2 ds 



^'11 I, will] l J -].. II t.,bfli! 2 I - / 

By lower semicontinuity of the norm and identification of the weak limit, we get the thesis. □ 



<i(||A 1 w(0)|| 2 + ||A 2 b(0)|| 2 )+ ^ (i(s),AMs))ds. 



5 Results for the second model 

In this section, we consider the following LES model for MBD, which is based on filtering only the 
velocity equation (and on the use of deconvolution operators): 

d t w + V • d (D Nl (w) ® D Nl (w)) - V • Gi (B <g> B) + - i/Aw = Gif, 
d t B + V • (B ® D Nl (w)) - V • (D Nl (w) ® B) - fiAB = 0, 
V-w = V-B = 0, 

w(0, x) = G lUo (x), B(0, x) - B (x), 

and we will work with periodic boundary conditions. A similar model in the case without decon- 
volution has been also studied in [9]. 

Here we take a 2 = 0, so that b = B and A 2 = G 2 = /, and 7V 2 = 0, so that Djv 2 B = IB = B. 
We set for simplicity 

a = ai>0, G = G U A = A u N = N%. 

The first aim of this section is to show the changes needed (w.r.t Thm. 3.1) to prove the existence 
of a unique solution to the system (5.1) for a given N G N, when we assume that the data are such 
that 

(5.2) u gH , BqGHo, and f G L 2 ([0, T] x T 3 ), 
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which naturally yields Giu £ H 2 , Gif £ L 2 ([0, T]; H 2 ). 

We start by defining the notion of what we call a "regular weak" solution to this system. 

Definition 5.1 ("Regular weak" solution). We say that the triple (w,B,q) is a "regular weak" 
solution to system (5.1) if and only if the three following items are satisfied: 
1 ) Regularity 

(5.3) weL 2 ([0,T];H 2 )nC([0,T];H 1 ), B £ L 2 ([0, T]; Hi) n C([0, T]; H ), 

(5.4) d t w G L 2 {[0,T];U ), d t B £ L 2 ([0, T}; H_i), 

(5.5) ggL^CT];^^)), 

Initial data 

(5.6) hm||w(f,-)-GiUo|| Hl =0, lim ||B(t, •) - B || Ho = 0, 

■2f\f\ Tl. IT 1 HIV. 1 31 



5; Weak Formulation : For all v,h £ L 2 ([0,T]; JJ 1 ^) 3 ), 
/0 Jt, 



(5.7) 



(5.8) 



f f <9 t wv- f f G!(%(w)®%(w)):Vv 

JO JT 3 JO JT 3 

[ [ Gi(B®B):Vv+/ / Vg-v 

Jo Jt 3 jo Jt 3 

v [ f Vw:Vv= / [ (Gif)-v, 

JO JT 3 JO JT 3 

8 t B h- [ [ (B®D Nl (w)):Vh 

j Jo JT 3 

/ / (D Nl (w)®B) : Vh + M / / VB : Vh 

JO JT 3 JO JT 3 



Remark 5.1. iJwe io i/ie certain symmetry in the equations, it turns out that B has the same 
regularity of D^w (not that ofw). 

All terms in the weak formulation are well-defined. Indeed, the only term to be checked (which 
is different from the previous section) is the bi- linear one involving B £ L 4 ([0, T]; L 3 (T3)) 3 and 
D N (w) £ L oo ([0,T];£ 6 (T 3 )) 3 . To this end, we observe that 

[ T [ (B®D N (w)) : Vh<G [ T \\B(t)\\ L3 \\D N (w)(t)\\ L6 \\Vh(t)\\ L2 

J0 JT 3 JO 

< Gr||B|| L 4([ 0T ] ;i3 )||£) A r(w)|| i oo([ 0:T ] ;i 6)||Vh|| L 2([ 0iT ] ;i 2). 

We have now the following theorem showing that system (5.1) is well-posed. 

Theorem 5.1. Assume that (5.2) holds, a > and N £ N are given. Then, Problem (5.1) has a 
unique regular weak solution satisfying the energy inequality 

|(||A*£>* (w)|| 2 + ||B|| 2 ) + H|VA^|(w)|| 2 + M ||VB|| 2 < G(||u ||, ||Bo||, iz-lf |U a([ o, n .H_ l ))- 

Proof. We use the same notation and tools from the previous section and the main result can be 
derived from the energy estimate. We just give some details on the estimates which are different 
from the previous case, since the reader can readily fill the missing details. We use Djv(w m ) in 
the first equation and B m in the second one as test functions to obtain 



' d (p 1 / 2 !)^ (wm)|| 2 + ||B m || 2 ) + v\\\7A 1 / 2 D 1 / 2 (vs m )\\ 2 +^||VB, 



|2 



2 dt 

= (A 1 /^V2 (Gfi/m)jA l/2^1/2 (wm)) ^ 
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Then, by using the same tools employed in the previous section, we have the following estimates. 



(5.9) 



Til 

Label 


Variable 


Bound 


(Jrder 


a) 


! 1 

A? £)?. (Wm ). B„, 


L°°([0,T] 


H )nL 2 ([0,T] 


Hi) 


0(1) 


b) 


£ ) J/ 2 (w m ) 


Z,°°([0,T] 


H )nL 2 ([0,T] 


Hi) 


0(1) 


c) 




L°°([0,T] 


H!)nL 2 ([o,r] 


H 2 ) 


0(a-i) 


d) 


w m 


L°°([0,T] 


H )n^([o,r] 


Hi) 


0(1) 


e) 


w m 


L°°([0,T] 


Hi) nL 2 ([o,T] 


H 2 ) 


0(0 


/) 


D N (w m ) 


L°°([0,T] 


H o )n^([0,T] 


Hi) 


0(1) 


9) 


D N (w m ) 


L°°([0,T] 


Hi)nL 2 ([0,T] 


H 2 ) 




h) 


<9 t w m 


L^([0,T];H ) 


o(«- x ). 


i) 


9*B m 


i 2 ([0,T];H-i) 


o(^T). 



The estimates (5.9-a)-(5.9-h) are the exact analogous of the corresponding ones from (3.11). 
What it remains to be proved is just (5.9-i). Let be given h £ L 2 ([0, T]; Hi); then 

(<9 t B m , h) = -m(VB to , Vh) + (B m ® i)jv(w m ),Vh) - (£>jv(w m ) 8> B m , Vh) . 

Hence we obtain, by the usual Sobolev and convex interpolation inequalities, 

|(fitB m ,h)| </x||VB m || ||Vh||+2||B m || £B ||D^w m || Z 8||Vh|| 

< ||VB m ||( M + Cp J vw„ 1 || 1 / 2 ||V^w m || 1 / 2 )||Vh||. 

Next, by employing estimates (5.9-a)-d)-f)-g), we get 



/ (d t B m ,h) dt 
Jo 

2 ([0,T];H 



< ||B m || L 2 ([0iT];Hl )(/i+ ||-Div(w m )||^ (f0T , H o) ||VDjv(w m )||^ aoT1H ) )||Vh|| i 2 ([0iT];i2) , 



.1/2 



'm;ilL~([0,T];Ho)l 

|Vh|| L 2([ 0T ] ;i 2). 



These estimates are enough to pass to the limit asm-4 +oo and to show that the limit (w, B) is 
a weak solution which satisfies 



(5.10) 



(5.11) 



I [ d t wv- [ I G{D N {w) ®D N (w)) : Vv+ / / G(B ® B) : Vv 

Jo Jt 3 Jo Jt 3 Jo Jt 3 

+ v [ f Vw:Vv= f / G(f ) • v, 
Jo Jt 3 Jo Jt 3 

d t B h- [ [ (B (g) Dn(w)) : Vh 
Jo Jt 3 

[ [ (-Djv(w) ® B) : Vh + ^i ( [ VB:Vh = 0. 

Jo Jt 3 Jo Jt 3 



T 

o Jr. 



The introduction of the pressure follows exactly as in the previous section, while the uniqueness 
needs some minor adjustments. Let in fact (wi,Bi) and (w2,B2) be two solutions corresponding 
to the same data (uo, Bo, f) and let us define as usual W := wi — w 2 and B := Bi B 2 . We will 
use ADn(W) and B as test functions in the equations satisfied by W and B, respectively. Observe 
that, by standard calculations, AD N {W) lives in L 2 {[0,T] x T 3 ) 3 , while B G L 2 ([0, T]; Hi). In 
order to justify the calculations — those for the velocity equation are analogous to the previous 
ones — first observe that, for any fixed order of deconvolution N, 
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since the duality is well-defined thanks to (5.9-a)-i). We formally write the distributional expression, 
keeping the time derivative, and we get the following equality (to be more precise, one should write 
directly the integral formula, after integration over [0,t], but the reader can easily fill the details): 

— (II^/^cmoII 3 + ll^ll 2 ) + Hiv^< 2 ( >v)|| 2 + m ||vb|| 2 

= -((D N (W) ■ V)D N (-wx), D N {W)) + ((23 • V)B X , D N (W)) 

- ((D N (W) ■ V)Bi, 23) + ((23 ■ V)£» JV (w 1 ) ) 23) 
= :h+h + h + h. 

Now, we need to estimate the four integrals in the right-hand side. The estimates are obtained by 
using the standard interpolation and Sobolev inequalities together with the properties of 2? at. We 
have: 

|/!| < e V \\VA^D^(W)f + Cs(N + 1)4SU ^° " VW1 " 4 \\A^D]i\w)f , 



\h\ < WBWuWVDxCW^WB^ 

<C||23|| 1 / 4 ||V23|| 3 / 4 ||V J D JV (>V)|| 1 / 4 ||A J D J v(>V)|| 3 / 4 ||B 1 || 

< c (N+l)^_ 1 1 g 1 1 1/4 1 1 VB | |3/4 a 1/4 1 1 v ^ 12 ( w) 1 1 l/4 a 3/4 1 1 A ^ 2 (W) || 3 / 4 ||B! || 

a 

^ (Ar + 1)V2 ||g|| l / 4 | | Vg | | 3/4|| A l/2^V2 (w)||1/4||VAl/2 ^l/2 (>v)||3 /4|| Bi || 



< eH|V23|| 2 + ev\\VA^D)l\wW + ffi^ff pil P 1/2 J# 2 (W)|| || Bl || 4 

< eH|V23|| 2 + ^IIV^/^OV)!! 2 + ^g±^!|| Bl || 4 (||23|| 2 + H^tfj/'OV)!! 2 ) 

IJal^llCivCWJIIz-llVBIIIlBxIl 

< C||V2? A r(W)H 1 / 2 ||A^ Ar (>V)|| 1 / 2 ||V23|| UBill 

< C ^+ 1 ) 1/2 ||AV^V 2 ( W) ||l/2||V^/2^ (n;) ||l/2 ||vg|| ||Bi|| 

< e/i ||V23|| 2 + Ce(7V 2 +1 V 1/2 < 2 (W)|| ||V^/^V2 (W) || || Bif 

< £/ ,|| V23|| 2 + eu\\WA^D]i\wW + ^ ^ Pi || 4 ||^ a jj/ a (W)|| a , 

|/ 4 | < ||B|| 2 4 ||V^(W)|| < C||S|| 1 / 2 ||V23|| 3 / 2 (7V + l)||Vwx|| 



::--HVfi||- + cr ' (jV ", + 1)4 ||Vw 1 n|g|- 



,2 , ^eV^jr^ IIY7,„ l|4|| H ||2 

M 3 

We then set e = 1/6 and, by collecting all the estimates, we finally obtain 



j^A^D^mf + ||23|| 2 ) + u\\VA l / 2 D)l 2 {W)\\ 2 + ,i||V23|| 2 < CM(\\A 1 / 2 d]( 2 (W)\\ 2 + ||23|| 2 ), 
where 

MMW+1 )Waup{e^,^^,i5lf;,M}. 

t > L ^ J /U d ii 6 / z v 6 l 2 or ii^voc ) 

An application of the Gronwall's lemma proves (for any fixed N) that P 1 / 2 D] V /2 (W)|| 2 + 1|23|| 2 = 0; 
hence, by using the properties of A and Dm exploited before, we finally get W = 23 = 0. □ 
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We can now pass to the problem of the convergence as N —¥ +00, proving the counterpart of 
Theorem 1.1. 

Theorem 5.2. From the sequence {(wjv, Bat, qN)}Nen, one can extract a sub-sequence (still de- 
noted {(wat, B^, gjv)}jVGN^ such that 

weakly in L 2 ([0, T]; H2), 
wat — »• w ^ weakly* in L°°([0,T]; Hi), 

strongly in £ p ([0, T]; Hi) Vpe[l,oc[, 

weakly in L 2 ([Q, T]; Hi), 
B \ — B < wea%* m i°°([0,T]; H ), 

strong in £ p ([0, T] x T 3 ) 3 Vp G [1, 10/3 [, 

g w — » g weakly m L 2 ([0, T}; H\T 3 )) n L 5/3 ([0, T]; VK 2 ' 5 / 3 (T 3 )), 

and sztc/i i/iai i/ie system 

<9 f w + V • G(Aw (g) Aw) - V • G(B <g> B) + Vo - i/Aw = f, 

V-w = VB=0, 

^ ' ' d t B + V • (B <g> Aw) - V • (Aw ® B) - ^Ab = 0, 

w(0, x) = Gu (x), B(0, x) = B (x) 
holds in distributional sense and the following energy inequality is satisfied: 

— iUwf + ||b|| 2 ) + H|VAw|| 2 + Ml|Vb|| 2 < (f, Aw). 
Proof. This result is based on the following estimates and from compactness results. 



(5.13) 



Label 


Variable 


Bound 


Order 


a) 


Wat 


L oo ([0,T];H )ni 2 ([0,T];H 1 ) 


0(1) 


b) 




L°°([0,T];Hi)ni 2 ([0,T];H 2 ) 


0(a- A ) 


c) 


D N (w N ), B N 


L°°([0,T];Ho)ni 2 ([0,T];Hi) 


0(1) 


d) 


qN 


L 2 ([0, Tli&fjs)) n £ 5/3 ([0, T]; W' 2 ' 5 / 3 (T 3 )) 


O^ 1 ) 


e) 




L a ([0,T];H ) 


0(0 


/) 


d t D N (w N ),d t B N 


L 4/3 ([0, T];H_i) 


0(1) 



The only new bound here is represented by the one for dtB^ from (5.13-f). In fact, by the usual 
interpolation inequalities, we get 

\(d t B N ,h)\ <m||VBjv|| [[Vh||+2[[B JV |U4||i?Ar(w i v)l| i 4llVhll 

< (^HVBjvII + C||B A ,|| 1 / 4 ||VB iV || 3 / 4 || J D 7 v(w 7 v)|| 1 / 4 ||V J D J v(w i v)|| 3 / 4 )||Vh||. 

Next, by employing estimate (5.13-c), we get 

\(d t B N ,h)\ < ( A1 ||VB N ||+C||VB J v|| 3/4 ||V^A r (wA r )|| 3/4 )||Vh||, 

and since both VBat, VDjv(w) g L 2 ([0, T]; L 2 (T 3 ) 9 ), we can show that 

i-T 

(d t B N ,h)dt < //||VB iV ]]i,2([ 0)T] ._ t ,2 ) ||Vh|| i 2 ([0)T] . I , 2) 

+ C||VBAr||^ 4 ([0 T];L 2 ) ||VL'Ar(w)Ar|| 3 : '2 4 ([0 T] . L2) \\Vh\\ L 4 {[0 T] . L2) , 

thus proving that d t B N G L 4/3 ([0, T]; H_i), independently of N. 

The limit N — > +00 can be studied as in the previous section. In addition to the same estimates 
proved before, from the bound on the time derivative of B we obtain that Bn — > B in L p ([0, T]; Ho), 
Vp G [1, 00 [, and reasoning as in (4.4) we get 

D N (w N ) (g> D N (w N ) — >Aw®Aw strongly in L p ([0,T] x T 3 ) 9 V>G[l,5/3[, 

Bat® Bat — > B ® B strongly in L p ([0, T] x T 3 ) 9 V>G[l,5/3[, 

D N {w N ) ®B N — >Aw®B strongly in L p ([0,T] x T 3 ) 9 V>e[l,5/3[. 

Finally, the proof of the energy inequality follows the same steps as before. □ 
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